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Abstract 



A number of recent experiments have showed that surfactants can modify 

the growth mode of an epitaxial film, suppressing islanding and promoting 

layer-by-layer growth. Here I introduce a set of coupled equations to describe 

the nonequilibrium roughening of an interface covered with a thin surfactant 

layer. The surfactant may drive the system into a novel phase, in which the 

surface roughness is negative, corresponding to a flat surface. 
PACS numbers: 68.55 -a, 68.35 Fx, 64.60 Ht 
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There is much theoretical interest in the statistical properties of nonequilibrium surfaces. 
The growing interfaces naturally evolve into self-affine structures; the surface morphology 
and the dynamics of roughening exhibit simple scaling behavior despite the complicated 
nature of the growth process In particular, much attention has focused on different 
models to describe thin-film growth by molecular-beam epitaxy (MBE) 

Under ideal MBE conditions the primary relaxation mechanism is surface diffusion, which 
conserves the mass of the film. Experimentally both lattice strain and surface free energy de- 
termine whether the film undergoes layer-by-layer growth, islanding, or layer-by-layer growth 
followed by islanding. In experiments involving growth of Ge on Si (100) surface layer-by- 
layer growth is limited to 3-4 monolayers (ML) due to the lattice mismatch between Si and 
Ge and is followed by formation of unstrained Ge islands. It was shown recently that is- 
landing in the Ge/Si system can be suppressed effectively by use of a surfactant monolayer, 
changing the growth mode from island growth to layer- by-layer growth |^. Suitable surfac- 
tants such as As and Sb strongly reduce the surface free energy of both Si and Ge surfaces 
and segregate at the surface during growth. A number of subsequent experiments showed 
that surfactants can change the surface morphology in a wide variety of systems: Sb altered 
the growth behavior of Ag on Ag(lll) |Q, antimony was found to change the structure of 
islands in Ge/Si growth P| and Te was used as surfactant to sustain layer- by-layer growth 
of InAs on GaAs(OOl) 0. 

Microscopically, arsenic placed on the surface of Si or Ge lowers the surface tension by 
eliminating the dangling bond states. But for ordinary surfactants reduced surface tension 
enhances roughness. The reduced adatom diffusion might change the growth mode, but so far 
most continuum models with or without surface diffusion predict rough interfaces P-|6|,p 



While the microscopic mechanism of As and Ge/Si interaction is much investigated [p!2 
understanding in the framework of continuum growth models is still lacking. 

In this paper I study the generic problem of nonequilibrium roughening of an interface 
covered by a thin surfactant layer (see Fig |T]). Building on experimental results and general 
symmetry principles, a set of nonequilibrium equations are proposed to describe the growth 



of an interface coupled to the fluctuations in the surfactant coverage. The analytic study 
of these equations indicates that the surfactant changes drastically the morphology of the 
interface in 2+1 dimensions. In particular, the coupled system supports the existence of a 
novel phase characterized by negative roughness exponent, which can be identified with a 
morphologically fiat surface. 

As mentioned above, under ideal MBE conditions, relaxation proceeds via surface diffu- 
sion. In contrast to the ideal MBE, there is experimental evidence that surfactant mediated 
growth of Ge on Si proceeds by highly local Ge incorporation with minimum surface diffu- 
sion 



131 . Ge atoms that adhere to the As-capped surface rapidly exchange sites with the As 
atoms and incorporate into subsurface sites. In the absence of surface diffusion, the growth 
equation may contain terms which violate mass conservation [Q. The simplest nonlinear 
growth equation with nonconserved dynamics was introduced by Kardar, Parisi, and Zhang 



(KPZ) 111: 



dth = z/V^/i + Xiyhf + r]. (1) 

Here h{x, t) is the height of the interface 'm d = d! + 1 dimensions. The first term on the 
right hand side describes relaxation of the surface by a surface tension v. The second term 
is the lowest order nonlinear term that can appear in the interface growth equation, and is 
related to lateral growth. ri{x,t) is a stochastic noise driving the growth; it can describe 
thermal and beam intensity fluctuations. Additional terms in ([^) will include the coupling 
to the surfactant fluctuations. 

An efficient surfactant must fulfill two criteria: it must be sufficiently mobile to avoid 
incorporation, and it must surface segregate. Careful experimental studies showed for the 
Ge/Si system that the bulk As concentration is less than 1%; thus the effect of As on 
growth is a surface phenomena |^. Neglecting the desorption of the surfactant atoms, 
the equation governing the surfactant kinetics obeys mass conservation. This leads to the 
continuity equation dtv = —V ■ j + where v{x,t) is the width of the surfactant layer 
rj' is a conserved uncorrelated noise which incorporates the random local fluctuations of the 



surfactant, and j is the particle-number current density. The simplest linear equation with 
conserved dynamics correctly incorporating the effect of surface diffusion is 

dtv = -KV\ + r}. (2) 

To account for the coupling between the growing surface and the surfactant it is necessary 
to introduce additional terms in Eq. (||) and @. There are two main criterias which restrict 
our choice: The coupling terms must satisfy the symmetry conditions characteristic of the 
interface and the obtained set of equations should be self-consistent, i.e. the resulting 
dynamics should not generate further nonlinear terms. In addition the coupling terms 
included in Eq. (Q) must obey the required mass conservation for the surfactant. 

The simplest set of equations that satisfy the above conditions is 

dth = vV^h + XiWhf + ^{Wvf + r/o (3a) 

dtV = -KV^v + iV\Vh) ■ iVv)] + 771, (3b) 

where the noise terms r^o and rji are assumed to be Gaussian distributed with zero mean 
and the following correlator: 

< r]^{x, t)ri,{x', t') >= V,5{x - x)5{t - t'). (4) 

Here 

Vo = Do (5) 

and 

Pi = -DiV^ + D2V\ (6) 
The D2 term is generated by Do and Di as will be shown below. 

The generic nonlinear term (Vf )^ in ( Pa| ) can be derived using symmetry principles. In 
(|bD the V^[(V/i) ■ (Vf )] term results from a current j = — V[(V/i) ■ (Vf )], and obeys mass 
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conservation. Geometrical interpretation of this term suggests that a positive 7 drives 
the surfactant to cover uniformly the irregularities of the surface, i.e. enhances the wetting 



properties [|T^. A negative 7 has the opposite effect, assigning a non-wetting character to 
the surfactant. Since in experiments there is no evidence of surfactant agglomeration (non- 
wetting character), but it is energetically favorable to terminate the Ge layer with As atoms, 
we assume that the surfactant wets the surface, thus 7 > 0. 

The quantity of main interest is the dynamic scaling of the fluctuations characterized by 
the width wl(t,L) = < [h{x,t) - h{t)]'^ > = L'^^'^ f {t / L""^) , where xo is the roughness 
exponent for the interface h{x,t), and the dynamic exponent zq describes the scaling of the 
relaxation times with the system size L; h{t) is the mean height of the interface at time t 
and the <> denotes ensemble and space average. The scaling function / has the properties 
f{u — s> 0) ~ M^^o/^^o and f{u — > 00) ~const. In a similar way one can define Xi^i-nd Zi to 
characterize the fluctuations in the surfactant coverage v{x,t). 

For /5 = 0, Eq. (^) reduces to the KPZ equation (|I|). For the physically relevant 
dimension, ci = 2 + 1, the scaling exponents are not known exactly. Extensive numerical 



simulations give xo = 0.385 ± 0.005 and zq ^ 1.6 |T^. Thus the interface is rough and the 
roughness increases with time as Wo(^) ~ t-^"'^^". 

For 7 = 0, Eq. (|3bD is the fourth order linear diffusion equation with conserved noise 



which can be solved exactly, resulting in zi = 4 and Xi = [19 



Thus, neglecting the coupling terms, Eq. ( pal) and ( pbD predict rather different values 
for z and the roughness exponents Xi- To see how the couplings change this behavior, I have 
investigated Eq. (^ using a dynamical renormalization-group (DRG) scheme. During the 
DRG calculations only one dynamic exponent z = zq = zi was used, valid if the system does 
not decouple. The fast modes are integrated out in the momentum shell e~'Ao <\k\< Aq, 
and the variables are rescaled as a; ^ e'x, t —>■ e^H, h —>■ e^°''h, and v e-^^'f . The 
calculations have been performed up to one-loop order. 

The diagrams contributing to A cancel each other, resulting in the scaling relation 
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^ + Xo = 2. (7) 

This relation is known to be the property of the KPZ equation and it is a consequence of 
Gahlean invariance (GI). Since the DRG conserves the GI, this scahng law is expected to 
remain valid to all orders of the perturbation theory. 

A second scaling relation can be obtained from the non-renormalization of the diffusion 
coefficient Di. 

z-2xi-d' -2 = 0. (8) 

The diagrams that contribute to Di have a prefactor proportional to k^, thus they are 
irrelevant (k is the wave vector in the Fourier space). They in fact contribute to D2, 
justifying its introduction in (^). 

These two scaling relations already indicate that the coupled interface/surfactant system 
is qualitatively different from the uncoupled one. For a planar interface {d' = 2) (|^,||) give 

Xo + 2X1 = -2, (9) 

which means that at least one of the exponents has to be negative. 

A third scaling relation unfortunately is not available, but insight can be obtained from 



numerical integration of the flow equations obtained from the DRG ||20|. The integration 
showed the existence of two main regimes: 

(i) In the flrst regime one or both of the coupling terms {(3, 7) scale to zero. In this case 
the two equations become completely (both coupling terms vanish) or partially (only one 
coupling term vanishes) decoupled, and the two equations might support different dynamic 
exponents z. The DRG scheme used is not reliable in this regime. 

(ii) The presence of a strong coupling flxed point is expected when both of the nonlinear 
terms diverge. The integration shows that this coupled phase exists only for z > 8/3 ||21|| . 
It is important to note that although there is no identifiable fixed point, in this phase the 
scaling relations (^ ^ are exact. According to (|^, for z > 8/3 the roughness exponent 
of the interface xo is negative. With a negative roughness exponent, every noise-created 



irregularity is smoothed out by the growth dynamics and the resulting surface becomes flat. 
Thus the coupling of the surfactant to the growing interface results in the suppression of the 
surface roughness. This corresponds exactly to the experimentally observed behavior, i.e. 
the addition of the surfactant suppresses islanding, resulting in a morphological transition 
from rough (without surfactant) to fiat (with surfactant) interface. 

The roughness exponent of the surfactant from (H) is negative if 2; < 4, while for z > 4 
it becomes positive. In the Ge/Si system, for example, the As has a saturation coverage 
of IML, which is independent of the system size and is governed only by the microscopic 
bonding of the As to the Ge dangling bonds. One expects no relevant fluctuations in the 
thickness of the coverage; this requires a negative roughness exponent for the surfactant and 
thus limits the dynamic exponent to values smaller than four. 

The DRG analysis fails to provide the exact value of the dynamic exponent z. As in 
the case of many other growth phenomena, simple discrete models might be very helpful to 



obtain its value |^ . Summarizing the results of the direct integration of the DRG equations, 
for z > 8/3 the existence of a strong coupling fixed point is observed, in which the interface 
roughness exponent is negative, corresponding to a fiat phase. There is no upper bound in 
z for the existence of this phase, but physical considerations suggest that 2; < 4, in order to 
allow the uniform surfactant coverage observed experimentally. 

It is important to note that introducing Eq. (^ I did not use directly the existence of 
the strain due to the lattice mismatch. Although an important problem [^], a continuum 



description of strain-induced roughening is still missing. The proposed model is expected 
to describe the coupled surfactant /interface system, but decoupling the surfactant does not 
necessary result in an equation describing heteroepitaxial islanding. Further studies are 
necessary to understand the microscopic (perhaps strain induced) origin of the nonlinear 
coupling terms. 

In Eq. d^) the desorption of the surfactant atoms is neglected by considering that (pE]) 
obeys mass conservation. Lifting the conservation law, (PB) should be replaced by a non- 



conservative equation. Such a system has been recently studied P4| , |25[| , and it was found 



that in most cases the couphng does not change the KPZ scahng exponents. Enhancement 
of the exponents is possible only when the couphng is one-way, i.e. one of the equations 
decoupled from the other one is acting as source of correlated noise. 

Further linear and/or nonlinear terms added to (|^) might influence the dynamics of the 
system. The goal here was to derive the simplest set of equations predicting the experimen- 
tally observed morphological phase transition; the study of other possible nonlinear terms 
and their relevance is left for future work. 

In conclusion, I have introduced a set of equations to describe the interaction of a growing 
surface with a surfactant. The main experimentally motivated requirements for @ were: (a) 
no surface diffusion of the newly landed adatoms; (b) conservative and diffusive surfactant 
dynamics, originating from neglecting incorporation and desorption of the surfactant during 
the growth process. The obtained equations indicate the existence of a coupled phase, in 
which two scaling relations between the three exponents are available. In this phase, the 
roughness exponent of the interface is negative, morphologically corresponding to a flat 
interface, as observed experimentally. 

Moreover, Eq. (|^) serve as a good starting point for future studies of an interface coupled 
to a local conservative field, a problem of major interest in the context of recent efforts to 
understand the general properties of nonequilibrium stochastic systems. 
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Sciences and the National Science Foundation. The Center for Polymer Studies is supported 
by National Science Foundation. 
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FIGURES 

FIG. 1. Schematic illustration of the studied surfactant/surface system. The figure represents 
a cross section of the two dimensional surface of heigh h{x,t) covered by a thin surfactant layer 
with thickness v{x,t). A newly arriving atom penetrates the surfactant and is deposited on the 
top of the growing interface h{x,t). 
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